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MAXIMUM  LIKELIHOOD  ANGLE  ESTIMATES  OF 
TWO  SIGNALS  USING  THREE  SQUINTED 
BEAMS  (SYMMETRIC  SOLUTION) 


INTRODUCTION 

The  problem  of  estimating  the  angle  of  arrival  of  closely  spaced  targets  has  received  considerable 
attention  (1-19].  Two  solution  methods,  the  maximum  likelihood  estimate  (MLE)  and  the  maximum 
entropy  method  (MEM)  have  been  used  extensively  for  many  problems.  For  array  antennas.  White  [1] 
calculated  the  mean-square  error  for  the  two-target  case  when  the  signal-to-noise  ratio  (SNR)  was  large, 
a  case  which  provides  small  errors.  This  calculation  was  made  for  both  the  case  in  which  the  plane  of 
symmetry  between  the  targets  was  known  and  the  case  in  which  it  was  not  known.  Lang  [2]  compared 
the  MLE  with  the  MEM  for  the  two-target  nonsymmetric  case.  He  found  that  the  MLE  was  only 
slightly  better  than  the  MEM;  at  most  the  difference  was  approximately  V2.  In  addition  to  comparing 
the  accuracy  of  the  MLE  and  MEM  for  the  two-target  case,  Trunk  et  al.  [3]  compared  the  MLE  accura¬ 
cy  to  the  Cramer-Rao  bound  for  unbiased  estimators.  They  showed  that  the  MLE  can  achieve  better 
accuracy  than  the  bound  for  the  lower  signal-to-noise  ratios  because  it  is  a  biased  estimator. 

Besides  the  accuracy  of  the  estimates,  a  number  of  implementations  have  been  studied.  White  [1] 
describes  a  double-null  tracking  system  which  approximates  the  MLE  estimate  using  an  array  antenna. 
He  also  discusses  a  tempered  double-null  tracker  using  multiple  beams.  Neither  of  these  implementa¬ 
tions  estimate  the  angle  of  arrival  on  a  single  pulse;  they  search  for  a  solution  through  a  tracking  loop. 
Howard  [4]  uses  the  MLE  method  to  obtain  the  angle  estimate  on  a  single  pulse  by  a  search  technique. 

Some  closed-form  solutions  for  the  angle  of  arrival  of  two  targets  under  special  conditions  have 
been  obtained.  Peebles  and  Berkowitz  [5]  obtained  a  noise-free  solution  for  two  targets  and  three 
squinted  beams  when  the  antenna  patterns  could  be  approximated  by  polynominals.  White  [1]  obtained 
a  fixed-beam  solution  under  the  special  conditions  of  symmetric  antenna  patterns  pointed  directly  at  the 
plane  of  symmetry  between  two  targets.  Recently  a  closed-form  MLE  using  three  subapertures  was 
obtained  for  two  targets  both  when  the  plane  of  symmetry  between  them  was  known  and  when  it  was 
unknown  [6].  The  accuracy  was  nearly  as  good  as  if  all  the  elements  of  an  array  had  been  used.  Since 
a  closed-form  MLE  using  three  subapertures  had  been  obtained,  we  addressed  the  analogous  problem 
of  obtaining  a  closed-form  MLE  using  three  squinted  beams.  We  first  reduced  the  mean  square  error 
cost  function  found  in  the  MLE  formulation  into  as  simple  a  form  as  possible  in  terms  of  the  un¬ 
known  angles  of  arrival.  We  then  found  the  angles  that  minimized  the  cost  and  evaluated  the  accuracy 
of  the  estimates.  We  investigated  only  the  case  in  which  the  plane  of  symmetry  separating  the  sources 
or  targets  is  known. 
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MLE  FORMULATION  FOR  THREE  SQUINTED  BEAMS  AND  TWO  TARGETS 

We  begin  by  defining  the  geometry  and  antenna  system  with  the  aid  of  Fig.  1 .  The  antenna  lies  in 
the  z-y  plane  with  the  principal  axes  of  the  antenna  defined  along  the  z  and  y  axes.  This  implies  that 
the  aperture  illumination  is  factorable  in  the  z  and  y  directions.  The  antenna  patterns  can  then  be  writ¬ 
ten  in  terms  of  the  product  of  two  individual  antenna  patterns  which  are  functions  of  the  direction 
cosines.  For  large  apertures  and  narrow  beams  pointed  near  the  bore  site  axis  x,  the  direction  cosine 
can  be  approximated  by  the  elevation  and  azimuth  angles  9  and  tf»,  respectively,  as  shown  in  Fig.  1. 
The  angles  of  arrival  of  the  radiation  from  source  1  and  source  2  are  defined  to  be  (@  j , «/*  o)  and  (92,4ioh 
respectively.  We  have  defined  the  azimuth  t/i0  of  the  two  sources  to  be  the  same  and,  consequently, 
there  is  no  cross-leveling  of  the  array  with  respect  to  a  line  connecting  the  sources. 


Fig.  I  —  Coordinate  system 


The  antenna  is  constructed  so  that  three  squinted  beams  are  formed.  The  beams  are  symmetric 
about  the  x-z  plane  and  are  squinted  in  the  elevation  direction  9 .  Using  representative  antenna  pattern 
shapes,  their  antenna  pattern  factors  in  the  9  direction  are  shown  in  Fig.  2. 


Fig,  2  —  Antenna  patterns  for  three  squinted  beams 
in  0  direction  for  a  fixed  angle  i b 
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The  three  beams  are  identical  and  the  separation  in  angle  of  adjacent  beam  centers  is  denoted  by 
0O.  The  pointing  of  the  three  squinted  beams  is  specified  by  the  pointing  angle  9C  of  the  center  of  the 
middle  beam.  The  plane  of  symmetry  between  the  two  sources  is  specified  by  angle  0,  —  (9 ,  -  02)/2. 
The  antenna  pattern  factors  for  the  three  beams  in  the  9  coordinate  are  given  in  terms  of  the  angle  of 
arrival  9  of  a  source  by  G(9  —  0t.  —  0O),  G(9  -  0C),  and  G(9  -  9C  +  0O),  respectively.  The  antenna 
pattern  factors  in  the  1 direction  are  given  by  //(«/»)  for  all  three  beams.  The  antenna  patterns  for  the 
three  beams  are  then  given  by 

F{9  - 9c-90,4i )  -  HU>)  <7(0  -0,  -0O). 

F(9  —  0(.,«fi)  =  HU>)  G(9  —  0,.), 
and 

F(9  -0C  +0O,«M  =  //(«!»)  G(0  -  9C  +0O). 

The  antenna  is  also  constructed  with  the  property  that  when  a  plane  wave  sweeps  across  the  aper¬ 
ture,  signals  will  be  produced  at  the  output  ports  of  the  beamforming  devices  that  are  all  in  phase. 
These  signals  after  conversion  to  baseband  are  denoted  by  the  complex  numbers  sh  s2,  and  s},  where 
the  real  and  the  imaginary  parts  are  the  time-coincident  samples  of  the  in-phase  and  the  quadrature 
channels  of  the  synchronous  detector,  respectively.  The  signals  sj,  s2,  and  s}  out  of  the  three  beam¬ 
forming  ports  can  be  expressed  as  the  summation  of  the  excitation  produced  by  the  plane-wave  signals 
sweeping  across  the  antenna  and  the  additive  thermal  noise.  For  two  sources  in  the  far  field,  the  re¬ 
ceived  signals  are  written  as 

S]  —  a,  F(.9 1  -  9C  -  0O,0)  +  a2  F(02  -  9C  -  0o,i/i)  +  n  |, 
s2  -  a\F(9 [  -  0f>)  +  a2  F(92  -  0r,tf»)  +  n2. 


and 


s3  -  a  |  F(9 ,  -  0C  +0O,0)  +  a2  F(92  -  9C  +0O,0)  +  n2, 


(2) 


where  a i  and  a2  are  the  strengths  of  the  sources  1  and  2,  respectively,  and  n i,  n2,  and  n}  are  the  ther¬ 
mal  noise  in  each  of  the  channels  associated  with  a  beamforming  device.  The  noise  is  assumed  to  be 
uncorrelated  and  Gaussian  distributed  with  a  mean  of  zero  and  a  standard  deviation  of  a .  Equations 
(2)  can  be  written  in  vector  form 


S-WA+N, 


(3) 


where 


«i 

«2 


and 


"II 

"12 

w  - 

"21 

"22 

"31 

"32 

The  coefficients  of  the  matrix  W  are  the  antenna  gains  of  each  of  the  beams  in  the  direction  of  each  of 
the  sources,  as  described  in  Eq.  (2),  and  are  given  by 


"II  " 

F{9 1 

-0,  -0fl,*), 

"12- 

F(92  ~  0,  —  0O.'1> ) 

"21  “ 

F(9 , 

-0,,<fr). 

"22- 

F{9  2  —  0,.,il'), 

"31  ” 

F(0, 

-  0,  +  0o.il'), 

"32- 

F(9  2  —  0,  +  0o,il') 

(4) 
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We  have  now  constructed  a  model  of  the  signals  received  by  an  antenna  using  three  squinted 
beams,  excited  by  two  distant  sources  at  different  angles.  We  would  now  like  to  estimate  the  unknown 
parameters  A  and  the  angles  of  arrival  0 ,  and  02  using  the  noisy  measurements  S.  The  estimation  tech¬ 
nique  we  will  use  is  the  maximum  likelihood  procedure. 


Since  N  is  a  Gaussian  random  vector  with  uncorrelated  but  equal  variance  components,  maximiz¬ 
ing  the  likelihood  function  is  equivalent  to  minimizing  the  square  error 

I  -  { S  -  WA)*  (S  -  WA),  (5) 

where  the  asterisk  represents  the  conjugate  transpose.  The  value  of  A  which  minimizes  L  is 

A  -  (WW)'1  W*S, 


which  expands  to 


A 


E~' 


(p i  -  ?i)  +  12W2  -  qi<?2> 
"21 


-(/>!  -  4i>  -  qi(P\Pi~  QiQj) 
"22 


Pi(Pi  ~  q7)  -  Q\(p\  -  q,) 
"21 


P\(P\  -  Q 1)  -  92p>2  -  ?2> 
"22 


-iPi-  «2>  -  Q\(P\Pi- 


"21 


(P2~  «2>  +  P\(PiP2~  «l«2> 


"22 


where 


E-  (p\-  q  1)2  +  (Pi  ~  Qi)2  +  (p\Pi  -  qm)2. 


"11  "12 

Pi  “  — -  « 1  “ 


q  2- 


"21  "22 

"Jl  "J2 

"21  ’  P 2  "22  ' 

Equation  (S)  for  the  mean-square  error  L  can  be  written  as 

1- 1  ii/i!2-  i  //•. 

1-1  i-i 


where 


for  / ' 


f  -  “  "n  <*1  “  "12  Pi 

1,  2,  and  3.  Substituting  the  values  a\  and  a2  found  in  (6)  into  (8)  we  obtain 
(Pi  -  qi)  (0>2  -  «2>*l  -  W2  -  *l«2>*2  +  (/>)  -  q  |)»3l 


fl  — 
/)- 


(P\  -  q |)2  +  (Pi  ~  qi )2  +  (P\Pi  -  «I«2>2 
-  (p\Pi  -  q\qi)  t(P2  -  qi)*\  -  (p\Pi  -  qi?2>»2  +  < />i  -  gi><?l 
(pi  -  «i)2  +  (pi  -  qi)2  +  (p\pi  -  W2)2 
(/>!  -  g|)  I0>2  -  ?2>*l  -  (P\Pl  -  »lff2>*2  +  0>l  ~  g|)^l 
(P\  -  q\)2  +  (Pi  -  «2>2  +  0>l/>2  “  ?l«2)2 


(6) 


(7) 

(8) 


(9) 
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Using  Eqs.  (9)  in  (5),  the  mean  square  error  becomes 

_  1 1  <P2  ~  92>S|  ~  (P\P2~  9|92>*2  +  (Pi  ~  <h)Sj\\2  (1{)) 

( P\  -  9()J  +  ( Pi  -  92>2  +  (PlP2  -  9i92>2 

The  coefficients  P\,  q i,  p2,  and  q2  in  (9)  are  complicated  nonlinear  functions  of  the  angles  of  arrival 
and  the  geometry.  Finding  the  parameters  which  minimize  the  error  would  be  difficult  in  most  cases. 
However,  a  special  case  is  examined  which  is  of  a  simpler  form  than  (10). 


SOLUTION  UNDER  SYMMETRY  CONDITIONS 


We  investigate  the  solution  for  the  angles  of  arrival  of  signals  when  the  angular  bisecting  plane 
between  them  is  known  and  both  emitters  are  at  the  same  azimuth  ■  The  peak  of  the  center  beam  in 
elevation  is  pointed  directly  at  the  bisecting  angle  between  the  two  sources.  In  the  three  squinted  beam 
cluster,  the  squint  all  appears  in  elevation.  For  these  conditions  9  f  -  0°,  9\  ”  —  02,  and  the  weights  be¬ 
come 


"ii  "  F(0 1  -  »o>). 

w12  -  F(-0 ,  -  Po.'f'). 

h>21  -  F(fix,^), 

w22-  F(-9X,$). 

•pji  ”  F<fi i  +P0.</'). 

Wj2-  F(-0 1  +P0,tfl). 

(11) 

Using  the  fact  that  the  antenna  patterns  are  factorable  under  certain  conditions  as  expressed  by  Eq.  (1), 
Eq.  (11)  can  be  rewritten  as 

w„-  //(*)  G(fix  -0O). 

w,2-  //(*)  G(-0x-0 0). 

w2X-  //(*)  G(P,), 

w22-  Hty)  CM i). 

w31-  //<*)  G(»,  +«0). 

h>j2«  //(d>)  G(—0 1 

(12) 

Using  the  definition  of  pu  qx,  p2,  and  q2  we  find 


*11  _ 

H(+)  G(9x-90) 

G(0,~0  0) 

Pi  “ 

"21 

H(+)  G(0X) 

G(9X)  ’ 

*31 

Hi*)  G(9X  +0O) 

G(0 1  +  0o) 

92" 

"21 

H(+)  G(9X)  G(0\)  ’ 

*12  _ 

H(*)  G(-fl,-fl0) 

G  (—0 1  -  0  o) 

9i  “ 

**'22 

//(*)  G(-0,) 

G(-0,)  ' 

"32  _ 

HU>)  (-0)  +90) 

G(—0 1  +  9  a) 

P  2" 

"22 

//(*)  GM,) 

G(-»,)  ’ 

Using  the  fact  that  the  antenna  patterns  for  the  beams  have  even  symmetry,  G(0)  -  G(-9),  we  find 
that  pi  -  p2  and  qx  -  q2  and  define 


and 


P  -  Pi  -  P2 


9  "  9i  "  92- 

Using  the  symmetry  conditions  in  the  expression  for  the  mean  square  error  L  in  Eq.  (10),  we  find 

,  IU|  -  Xs2  +  s3||J 
2  +  X2 


(13) 
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,  a  _ —  — 


where 


X  -  P  +  Q 


G(0|  -  0O)  +  G[-(0,  +  0O)1 
G(0,) 


To  minimize  the  mean  square  error  L  with  respect  to  the  angle  of  arrival  0,  in  (13)  is  equivalent 
to  minimizing  L  with  respect  to  X  and  then  relating  X  to  0  \  through  a  table  look  up.  The  minimum 
value  of  L  is  found  by  finding  its  stationary  points  with  respect  to  the  real  variable  X  in  terms  of  the 
complex  numbers  st,  s2,  and  s3.  The  stationary  points  are  found  by  differentiating  the  mean  square  er¬ 
ror  L  with  respect  to  X,  setting  the  result  equal  to  zero,  and  solving  for  X.  The  result  is 

*  =  k  ~  hi  ±  /  (v  -  M  7 2,'  (i4) 

a  +  /3  V  a  +  /3 

where 

a  =  s,  *  s2  +  s2*  s,. 

/3  -  s2*  s3  +  s,*  s2, 
y  —  s2*  s2,  and 

TJ  —  J|  J|  *  +  s3  Jj*  +  S |  Jj*  +  s}  S |*. 


Because  of  the  geometry,  0 ,  cannot  be  less  than  zero  degrees,  and  consequently 

G(-0  p)  +  G(—0  0) 

X  "  C(0) 


or 


X  > 


2G(90) 
G( 0)  ’ 


where  G(0n )  ”  G(-00).  If  X  is  less  than  this  zero-angle  value  due  to  noise,  we  limit  X  to 

„  2G(9q) 

G( 0)  ’ 

The  solution  is  that  value  of  X  given  by  (14)  or  (15)  which  makes  L  in  (13)  a  minimum. 


(15) 


Not  only  can  the  angle  of  arrival  0,  be  estimated,  but  the  estimates  of  the  complex  amplitudes  can 
be  made  using  (6).  For  the  special  case  under  consideration,  a |  and  o2  are  estimated  in  terms  of  pand 
<7  by 

1  [1  +  p(p  +  ?)]  S|  +  (p  -  q)s2  -  [1  +  q(p  +  <?)1  s3 

0|  w,2  ( p  -  q)  [(p  +  q)2  +  2] 

and 

1  -  [1  +  q(p  -1-  q)  1  S]  +  (p  -  q)s2  +  [1  +  pip  +  <7)]  s2 

w,2  ( p  -  q)  (( p  +  <7 ) 2  +  2] 


To  compare  the  results  to  previous  work,  we  define  a  reflection  coefficient  by 

a2  -  [1  +  q  (p  +  (7>1  +  (p  -  q)s2  +  [1  -f -  p  (p  +  q)  s , 

P  <?,  (1  +  p  (p  +  <7)]S|  +  (p  -  q)s2  -  [1  +  q  (p  +  <7))$, 
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If  the  angle  of  arrival  and  stored  antenna  patterns  are  known,  p  and  q  can  be  found.  Consequently,  the 
reflection  coefficient  can  be  estimated.  For  the  latter  purposes  we  define  the  magnitude  Ip  I  and  phase 
</>  of  the  reflection  coefficient  p. 


For  a  system  to  be  put  in  operation,  a  calibration  curve  of  X  vs  the  angle  of  arrival  9 !  would  be 
found  experimentally  and  stored  for  reference  in  a  table.  For  our  purposes  here,  we  wish  to  illustrate 
the  nature  of  this  curve  and  use  it  later  to  indicate  the  system’s  expected  performance.  The  normalized 
antenna  patterns  chosen  are 


G(9  ±  0O) 


sin  nd(9  ±  9 0)/X 

it d(9  ±  90)/\  i t2  —  [ird(9  ±  90)/\]2 ' 


(17) 


where  d  is  the  aperture  width  and  X  is  the  wavelength.  The  beamwidth  Be  for  this  antenna  pattern  is 

Be  -  82.5 k/d 


in  degrees.  For  an  antenna  aperture  width  of  2.5  wavelengths  and  a  squint  angle  of  0.8  beamwidth,  the 
calibration  curve  of  the  parameter  X  vs  angle  of  arrival  0|  using  the  postulated  antenna  pattern  (17)  is 
shown  in  Fig.  3.  We  find  that  the  curve  is  very  flat  for  small  angles,  which  indicates  that  large  signal- 
to-noise  ratios  will  be  required  to  obtain  descent-angle  estimates  for  closely  spaced  targets.  We  next 
evaluate  quantitatively  the  performance  of  the  angle-measurement  procedure. 


Fig.  3  —  Calibration  curve 
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ANGLE  ESTIMATE  PERFORMANCE 


The  performance  of  the  MLE  angle  estimates  has  been  evaluated  using  a  simulation  technique. 
To  compare  the  results  to  previous  work,  the  signal-to-noise  ratio  is  defined  to  be 


N 


10  log 


kl 

2/t  ‘ 


This  is  the  signal-to-noise  ratio  of  one  source  (the  strongest)  at  the  center  of  the  center  beam.  Using 
this  signal-to-noise  ratio,  noisy  signals  are  generated  using  (3),  and  the  angle  of  arrival  9,  is  estimated. 
This  process  is  repeated  400  times  and  the  root  mean  square  (RMS)  error  in  the  estimates  of  0U  |p|, 
and  <f>  is  computed. 


For  a  30-dB  signal-to-noise  ratio,  a  magnitude  of  the  reflection  coefficient  of  0.9  and  three  values 
of®,,  the  RMS  error  in  estimating  9 ,  is  plotted  vs  the  relative  phase  shift  <j>  between  the  sources  in  Fig. 
4.  We  find  that  the  error  is  smallest  when  the  effective  signal-to-noise  ratio  is  largest  (i.e.,  signals  a, 
and  a2  are  in  phase)  and  largest  when  the  signals  are  out  of  phase.  Also,  Fig.  4  shows  that  the  closer 
the  signal  spacing,  the  more  difficult  it  is  to  measure  the  angle  of  arrival.  The  RMS  error  in  estimating 
the  angle  of  arrival  is  also  shown  in  Fig.  5  using  the  same  conditions  as  used  in  Fig.  4,  except  that  | p  | 
=  0.2.  The  effective  S/N  is  more  constant  with  relative  phase  between  the  sources,  and  so  the  curves 
are  flatter.  Again,  Fig.  5  shows  that  it  is  easier  to  estimate  the  angle  of  widely  separated  signals. 


(ft 


PHASE  DIFFERENCE  BETWEEN  a,  AND  a2  f  (dtgrMa)  PHASE  DIFFERENCE  BETWEEN  a,  AND  a2  A  (<U«rMt) 


Fig.  4  —  RMS  error  in  angle  of  arrival  for 
Ip  I  -  0.9  and  S/N  -  30  dB 


Fig.  5  —  RMS  error  in  angle  of  arrival  for 
Ip  I  -  0.2  and  S/N  -  30  dB 


Besides  estimating  the  angle  of  arrivals,  we  estimated  the  relative  amplitude  and  phase  between 
the  two  signals.  Using  a  30-dB  signal-to-noise  ratio  and  a  reflection  coefficient  magnitude  of  0.9,  the 
RMS  error  for  |p  I  and  <t>  is  shown  in  Figs.  6  and  7,  respectively,  for  three  angles  of  arrival.  Again  we 
find  that  the  estimates  are  better  for  the  more  widely  separated  sources. 
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Fig.  6  —  RMS  error  in  the  magnitude  of  the  reflection  Fig.  7  —  RMS  error  in  the  phase  of  the  reflection  coefficient 

coefficient  for  Ip  I  =  0.9  and  S/N  =  30  dB  for  Ip  I  =  0.9  and  S/N  =  30  dB 


We  next  compare  the  performance  of  the  MLE  solution  using  three  squinted  beams  to  the 
Cramer-Rao  bound  results  found  in  Ref.  3  and  to  the  MLE  solution  using  three  subapertures  16] .  This 
comparison  is  shown  in  Fig.  8,  where  a  30-dB  signal-to-noise  ratio.  Ip  I  =  0.9,  and  two  different  angles 
of  arrival  are  used.  We  find  the  errors  in  the  MLE  solutions  using  either  subaperture  or  squinted 
beams  are  not  very  different.  Sometimes  better  results  than  predicted  by  the  Cramer-Rao  bound  for 
unbiased  estimators  can  be  achieved  with  biased  estimators  than  with  unbiased  ones.  Although  not 
shown,  the  estimates  for  Ip  I  and  are  about  the  same  for  the  MLE  estimates  using  subapertures  [61  or 
using  squinted  beams  as  described  in  this  paper. 

SUMMARY 

The  estimation  of  the  angle  of  arrival  of  two  closely  spaced  sources  is  studied  using  an  antenna 
constructed  to  form  three  identical  squinted  beams.  A  simple  closed-form  solution  for  the  angle  of  ar¬ 
rival  is  found  under  certain  restrictions.  The  restrictions  are  that  the  center  of  the  middle  beam  is 
directed  toward  the  bisecting  line  between  the  targets  and  that  the  principal  axes  of  the  antenna  are 
parallel  to  and  perpendicular  to  the  line  connecting  the  targets.  This  condition  can  be  obtained  with  an 
antenna  viewing  an  emission  or  reflection  over  a  smooth  reflecting  sea. 

The  accuracy  of  the  estimates  of  the  angle  of  arrival  and  the  reflection  coefficient  was  studied  us¬ 
ing  simulation  procedures.  The  best  angle  estimates,  especially  for  reflection  coefficients  near  I,  weie 
obtained  when  the  signals  were  in  phase  and  the  worst  when  the  signals  were  out  of  phase.  This  can  be 
explained  by  the  high  and  low  effective  signal-to-noise  ratios  under  these  conditions.  The  results 
showed  the  estimates  became  better  as  the  target  separation  increased.  Finally,  the  MLE  estimates  for 
three  squinted  beams  were  comparable  to  the  MLE  estimates  for  three  subapertures;  only  minor 
differences  were  noted. 
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PHASE  DIFFERENCE  BETWEEN  a,  AND  at  +  (d*gr«*«) 

Fig.  8  —  Comparison  of  RMS  angle  of  arrival  errors 
for  Ip  I  -  0.9  and  30  dB  S/N 
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